The role of environmental correlations in the non-Markovian dynamics of a spin system 
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We put forward a framework to study the dynamics of a chain of interacting quantum particles affected by 
individual or collective multi-mode environment, focussing on the role played by the environmental quantum 
correlations over the evolution of the chain. The presence of entanglement in the state of the environmental 
system magnifies the non-Markovian nature of the chain's dynamics, giving rise to structures in figures of 
merit such as entanglement and purity that are not observed under a separable multi-mode environment. Our 
analysis can be relevant to problems tackling the open-system dynamics of biological complexes of strong 
current interest. 
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The investigation on the interaction between a spin-like 
two-level system (TLS) and an environment embodied by an 
ensemble of harmonic oscillators, the so-called spin-boson 
model 0]], has gained renewed momentum due to the cur- 
rent strong interest in the quantum-mechanical description of 
excitation-transfer processes in biological organisms or the 
dynamics of protein-solvent systems in photosynthetic com- 
plexes |0-@1 ■ The interest in such a celebrated model is 
twofold, in this context. On one side, it is employed to un- 
derstand the active role of an environment in assisting trans- 
port across such structures ]4[]. On the other hand, the model 
is extremely useful to characterize the extent to which non- 
classical behaviors survive in systems affected by environ- 
mental worlds that are strictly adherent to the paradigm of 
being wet and warm. In this respect, we aim at gathering a 
description of the TLS -environment evolution that is as com- 
plete as possible and able to account for subtle yet potentially 
key features of such dynamics emerging from the finiteness 
the environment itself, its capability to retain information on 
the system it is coupled to (and kick it back at later times) and 
its inherent structure. For, instance, Rebentrost et al. have 
shown in Ref. |01 that environmental memory-keeping effects 
are crucial for preserving the coherent behavior of excitonic 
energy transfer processes in biological complexes. 

Stepping aside from interesting yet still controversial bio- 
physical scenarios, the general model addressed above plays 
an important role in various contexts of applied quantum 
physics, from cavity quantum electrodynamics to its solid- 
state analogues involving defects in diamond coupled to semi- 
conducting microcavities HI] or superconducting devices and 
planar strip-line resonators J^]. In the latter case, very in- 
teresting regimes are quickly becoming available where the 
coupling strength of the coherent part of an evolution can be 
made comparable to the one responsible for the open-system 
dynamics due to the presence of the environment iflfJl . 

All such considerations and physical configurations de- 
mand the development of a framework apt to encompass 



and harness features of non-Markovianity in various coupling 
regimes. Formally, one is required to retain, to some extent, 
the degrees of freedom of the environment so as to allow for its 
possibly-non-trivial structure to influence the dynamics of the 
system of interest. This can be computationally very demand- 
ing and while sophisticated techniques have been put forward 
to determine the quantum statistical properties of interacting 
quantum many-body systems, they are either hardly suitable 
for non-linear configurations or face some difficulties in track- 
ing the temporal dynamics of the system itself. 

In this paper we contribute to the investigation on the in- 
teraction between a system of interacting TLS's coupled to a 
bosonic environment by studying the effect that environmen- 
tal entanglement has on the dynamics of a chain of interacting 
TLS's. We build up a mathematical apparatus that, inspired by 
and significantly extending the approach in Ref. fllll . solves 
exactly the local interaction between a single TLS and its 
multi-mode environment, so as to build an eigenbasis that 
is then used in order to tackle the inter-particle coupling by 
keeping track of the environmental properties. This helps us 
in incorporating the full quantum features of the environmen- 
tal state, therefore enabling a quantitative study of the role of 
entanglement in the evolution of the chain under investigation 
and its relation to the non-Markovian nature of the dynamics. 
We illustrate our techniques by addressing the cases of a dimer 
and a trimer, which are the smallest non-trivial configurations 
allowing us to study the most crucial attributes of such a com- 
plex problem. Our investigation opens up the possibility to 
study full-scale networks of coupled TLS's in an open-system 
fashion, covering the full range of coupling strengths with the 
environments, from the weak to the strong one. Although in- 
evitably leaving a few questions open, we believe such an en- 
deavor to be helpful in the contexts described above. 

The remainder of this paper is organized as follows: in 
Sec. U we describe the method put forward in order to ana- 
lyze the dynamics of an iV-site chain interacting with indi- 
vidual multi-mode bosonic environments without resorting to 



2 



the weak-coupling regime and thus allowing for a more com- 
plete overview of the effects of non-Markovianity on the evo- 
lution of the chain. Sec. |IT]is devoted to the explicit study on 
a dimer of interacting particles affected by entangled environ- 
ments. The key features in the dynamics of such configura- 
tion are compared to what is achieved under the assumption 
of independent environments so as to show that entanglement 
appear to enhance the non-Markovian aspects of the time evo- 
lution. Sec. Hill addresses the case of a common environment 
affecting the sites of our dimer, showing that the manifesta- 
tions mentioned above are insensitive to the non-local nature 
of the environment, but rather depend on the properties of its 
state such as the entanglement-sharing structure, as studied in 
Sec. [lV]for the case of a trimer. Finally, in Sec. [V] we draw 
our conclusions and briefly discuss the questions opened by 
our investigation. 



I. THE MODEL AND THE ANALYTICAL METHOD TO 
ITS APPROACH 

The model under consideration is described by an Hamilto- 
nian of the form 

H to t = Hs + Hss + Hb + Ssb- (1) 

Here Hg = —fi^2jLi e j&j gives the free dynamics of 
the elements of the chain, Cj being the energy split be- 
tween the levels of each TLS, while {<rj, fij, ct| } are 
the Pauli matrices for particle j=l,..,N. Analogously, 

Hb= Z)j=i Sfeli f^j,kaj fc Oj,fc is the free energy of the en- 
vironmental harmonic oscillators. The pedex k labels the rij 
independent modes interacting with the TLS at site j. Each 
of such modes has frequency uij^ and is described by the 
bosonic annihilation [creation] operator a,j ; k [fit fe ]. On the 
other hand, Hss models the inter-particle coupling within the 
chain, while Hsb describes the TLS -environment interaction. 
We take these two terms as 

N nj 

Hsb =Y, & j ® Yj 9j,k( & j,k + 

j=1 (2) 

JV-l 

Hss = - 2h ]T J jij+1 (ap] +1 + ap] +1 ) 

n=j 

with Jj.j+i being the nearest-neighbor coupling strength be- 
tween the particles occupying sites j and j + 1 on the chain 
and gj t h the analogous parameter associated with the interac- 
tion of the j th TLS and its fc th environmental mode. 

The first of Eq. generates a displacement of the environ- 
mental oscillators, conditioned on the state of the correspond- 
ing TLS. Our approach to the description of the dynamics is to 
solve the problems embodied by each the subsystem made out 
of a single TLS and the ensemble of local modes it interacts 
with. We thus consider fi = fl s + fi B + H S b and look for 
eigenstates of the form \a)j <8> \<j>i)j with a = ± and {|±) •} 
denoting the eigenbasis of <r|. States \(f>i)j are the eigenstates 



of the total Hamiltonian of an ensemble of oscillators (per as- 
signed state of the corresponding TLS), that is 

^2[±gj,k( a ik+ a j,k) + ^filfcaj,fc|0±>j = (-B 3 ± tj)\<l>±)j, 

k=l 

(3) 

where is the total energy of the j th subsystem. A detailed 
calculation shows that such energy depends on the number of 
excitations n^. in the oscillators coupled to the TLS at the j th 
site according to the expression 

K = ±( j+Y. to,* (nj,k - . W 
fe=i 

where atjk = gj,k/fiUj,k- The corresponding eigenstates 
can be written as \ip^ )j = \±)j®\a±)j=\±)j®^ 1 \a± h ) j , 
where rij = {ni,ri2, ..,n nj } is the vector specified by the 
number of excitations populating the rij oscillators interact- 
ing with the j th TLS and we have introduced the excitation- 
added coherent states |aj fc ) = [(a^ ± a :( - ; fc)™ fc /v / nfeT|| ± aj,fc)- 
Although it is straightforward to check that (a+ |a")=<5 m „ 
((a m |a'l)=(5 m „)and the set of states {|+,o^)} ({|-,a™}}) 
form an orthonormal basis, states |a!L) and |a") are not 
mutually orthogonal. A treatment with an arbitrary number 
of modes per site, performed along the lines of the discus- 
sion above, is certainly possible although computationally de- 
manding. Therefore, from now on, we will focus our attention 
to the case of a single mode per site and drop the mode-label. 
In the oscillator basis discussed above, the Hamiltonian (Q]) 
takes the form 

N 

H tot =]Th + Hss, (5) 

i=i 

each Hq being now diagonal. The remaining term Hss, 
which is defined in the second line of Eq. (0, only contains 
the degrees of freedom of the TLS chain. In the following, 
we consider a homogeneous chain with Jjj+i = J, Vj = 
1,..,N, By introducing the raising and lowering operators 
frj = (erj ± icrJ)/2, it is convenient to rewrite it as 

N-l 

Hss = -%J J2 + <V^ + +i)- < 6 > 

n=l 

In the interaction picture defined with respect to the unper- 
turbed Hamiltonian flo= J2j H i we get 

(N \ N-l I N \ t 

n &t e « °n+i + h.c.) n 
3 = 1 J n=l \fe=l / 

(7) 

withe± = Er= ^ st l±,«±)(±,«±l- 
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After elaborating on this expression, we obtain 

N-l 

Hi(t) = -hj (a+a- +1 e+e- +1 e-te+| 1+ ^. c .)- (8) 

n=l 

By expanding the O^'s in a coherent-state basis, we finally 
arrive at the effective expression for the intra-chain Hamilto- 
nian 



N-l 



n+l 



r~r. 



n+lJj 



(9) 



where = ct^ ® Dj[±a(i))] with D(£) the displacement 
operator of amplitude £ g C 02] and ay (i) = 2a J -(l-e* a ''*). 
While the analysis above explicitly consider individual envi- 
ronments affecting the elements of the chain, one can also 
study the complementary situation where all the TLS's col- 
lectively interact with a single (and common) ensemble of os- 
cillators. In this case, [Htot> Y^j 0«] = O> which in turn implies 
that Hss can be easily diagonalized. We tackle an explicit ex- 
ample of this situation in Sec. [HI] 

In order to obtain information on the full dynamics of the 
TLS subsystem, we write the density matrix p(t) of the chain- 
oscillators system in the excitation-added displaced basis. For 
all our simulations, we assume the factorized initial state 



P(0) = /Ochain(O) <g> Pho(O) 



(10) 



with /9chain(t) [Pho(t)] describing the state of the chain (oscil- 
lator environment) at time t. In what follows we consider two 
significant instances of Pho(O): the tensor product of individ- 
ual thermal states and the case of entangled oscillators pre- 
pared in iV-mode squeezed states 02] ■ As it will be discussed 
later on, such three paradigmatic cases allow us to well iden- 
tify the effect of the entanglement among the environments on 
the dynamics of the TLS's. 

Let us start addressing the case of a thermal state 



N 



t — 1 /— n J > 



3=1 1=0 



where rij={eP tiJJJS — is the average number of quanta in 
the oscillator at temperature T (we have taken the inverse tem- 
perature /3 = 1/fchT with kf, the Boltzmann constant) and 
|Z) is the state with I excitations. In the excitation-added dis- 
placed basis of the oscillators, we can recast p± into the form 



A' 



Pth z 



E 
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= 1 n,m,l=0 y ■" 



(12) 

where the + or — sign for the j th oscillator has to be chosen 
according to the state of the corresponding element in p c h a in- 
The scalar products (a±\l) needed for the evaluation of the 
expression above are given in terms of Charlier polynomials 



C„(m;a 2 )as01[il 



a n+m e -a /2 

(a"|m) = — C n (m;a 2 ), 

<\/m' 



'n\ym\ 

a n+m e -a 2 /2 

(a+|m) = C m (n;a ), 

n\\ m' 



(13) 



where we have assumed a g K. In the remainder of the 
manuscript we take Hj — n, Vj. We will also consider an 
TV-mode a pure squeezed- vacuum state for the harmonic oscil- 
lators \sq N ) =(l-A 2 )- 1 /2 J2 t \ l (S)|Li with A = tanhr 
and r the squeezing parameter. This state can be expressed 
in the chosen oscillator basis following the very same lines 
highlighted above and using again the Charlier polynomials. 
It is worth emphasizing that these two environmental states 
are locally indistinguishable as a squeezed vacuum is locally 
equivalent to a thermal state. 



n. CHAIN DYNAMICS UNDER INDEPENDENT 
ENVIRONMENTS 



We are now in a position to start addressing the effective dy- 
namics of the TLS chain under the influences of the various in- 
stances of environment introduced in the previous Section. In 
order to strip down our discussion from unnecessary compli- 
cations, we study here the smallest configuration of the model 
at hand, i.e. a two-TLS chain (a dimer) affected by two indi- 
vidual harmonic oscillators, one per site. The two TLS will 
be considered as identical. While embodying an interesting 
enough situation [2], this case allows for an agile description 
of the physics involved. 

In the dimer basis, the interaction hamiltonian in Eq. (O 
takes the form 



Hi(t) 



/ 


^(01)^2(02) 

V 



\ 

Di(oi)Dj(a 2 ) 

J 



(14) 



In order to simplify the notation, we take identically coupled 
TLS-oscillator subsystems, so that ctj = a (j=l,2) (this 
assumption, as well as the one of identical TLS's, does not 
limit the generality of our results and can be relaxed with only 
mild complications). The corresponding time-evolution op- 
erator U(t) is then calculated as a time-ordered Dyson series 

U(t)=T, m Un(t) With 



u m {t)-- 



(-i)> 



t pt 

dti / dt2- 
Jo 



d^fHi^Hj^.-Hi^), 
(15) 

where 7 is the time-ordering operator. The anti-diagonal form 
of Hi(t) helps considerably in working out a manageable 
expression or the time-evolution operator. In fact, the even 
(odd) terms in the expansion of U(t) involve diagonal (anti- 
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FIG. 1: (color online): Coherences of the dimer's density matrix 
for a chain affected by an environment prepared in a thermal and a 
squeezed- vacuum state [panel (a) and (b), respectively]. We have 
taken identical parameters for the two TLS-oscillator subsystems 
with hJ = 0.04eV, hw = 0.02eV, g = 0.02eV at room tempera- 
ture (T = 293 K). In this case a = 1. As explained in the text, the 
squeezing factor for the results in panel (b) has been taken so that the 
reduced single-oscillator state has the same effective temperature as 
in the thermal state case [i.e. panel (a)]. Quantitatively, r = 0.816. 
A coherence revival at t — n/oj is present for a squeezed- vacuum 
environment: the presence of quantum correlations between the en- 
vironmental modes fundamentally affects the dynamics of the TLS's. 



diagonal) operators, in the dimer basis. In details, we get 



2m 



5>n©+-(o(*3fc-i))^_(a(tafc)) 



fe=i 



fe=i 



2<- + l®n^-+( Q (*2fc-l))^-+(c^2fc)), 
fc=l 

m 

I fti(**)=| + -)i2<- + \®l[V + -(a(t 2k ))vl + (a(t 2k+1 )) 



such two different environmental preparations. Fig. [2] shows 
the purity P(t) = Ti(pf 2 ) of the chain's density matrix p 12 
in both the thermal and squeezed case. We plot such a fig- 
ure of merit against the intra-site coupling strength J and the 
re-scaled interaction time u>t. Quite reasonably, at values of 
J that are comparable with the TLS-oscillator coupling, the 
differences between the thermal and squeezed-vacuum cases 
becomes very small. We can understand this by considering 
that, for decreasing values of J, the non-local interaction be- 
tween the sites of the chain becomes less relevant than the 
local TLS-oscillator coupling. This implies that each TLS is 
strongly tied to its associated mode so that the local features 
of the oscillator environment become more preponderant than 
the non-local ones and the differences with a configuration 
of independent environmental states fades away. On the con- 
trary, the inter-site tunnelling arising from a large value of 
J makes the inter-mode entanglement relevant. Such differ- 
ences manifest themselves in a rather distinct behavior of the 
properties of the TLS state. 

As a special but very interesting case, we consider the 
chain to be initialized in the maximally entangled state 

(l/v2)(| ) + | + +))i2 an d wonder about the behavior 

that entanglement has under the interaction model we are 
studying. Such an initial state is interesting on its own as, 
remarkably, it is relatively straightforward to find closed ana- 
lytical expressions for the elements of the reduced TLS den- 
sity matrix pi 2 . While the populations of such matrix remain 
constant, the interesting element turns out to be the coherence 
(+ + |/3i2 (t) | ), for which we find 



2m+l 



(+ + \pi2(t)\--)=-e 



1 4»£t+8a 2 [l-co8(ait)]t 



(17) 



k=l 



fe=0 



for an oscillator environment prepared in the tensor product of 
two individual thermal states, while 



k=0 

where V + -(a(t k ))=Di(+a{t k ))D\(-a{t k )) [with an anal- 
ogous definition for T> ^(a(t k ))]. 

By initializing the first (second) TLS in the coherent super- 

with \a\ 2 + |6| 2 =1 



in 



=o - i + 



. 1 4»£t+8a 

\pi 2 (t)\--)=-e 



2 H"^S, + 2A sq coa(u)t) 



(18) 



in r, 



0, the state 



|-} 2 ) and considering that [H to t> X^=i ^j) 
\+, +)i2 will stays unpopulated, while the probability to pop- 
ulate | — , — ) 12 remains constant. In Fig. [T] we plot the co- 
herences of the density matrix of the two TLS's affected by 
boson environments prepared in either uncorrected thermal 
states or a squeezed-vacuum one. In order to perform a faith- 
ful comparison between such two cases, we have taken the 
squeezing parameter according to tann r = e 2 with the 
value of j3 provided by the thermal-state environment. This 
ensures that the reduced single-oscillator state obtained from 
the squeezed-vacuum environment is a thermal state having 
the same temperature as the one considered in the genuinely 
thermal instance. In this way, we can nicely isolate the effect 
of the entanglement shared by the environmental oscillators. 

Clearly, the revival of coherences at t = ir/uj is only due 
to the presence of correlations in the environment and an in- 
dication of a significant deviation of the state properties in 




FIG. 2: (Color online): Purity of dimer's density matrix affected by 
individual environments prepared in a thermal state (red surface) and 
a squeezed-vacuum one (green surface, staying above the previous 
one). P(t) is plotted against the chain inter-site coupling energy hJ 
for hu = 0.04eV, g = 0.02eV and T = 293K, which corresponds 
to a — 1/2. We have taken a squeezing factor r = 0.489 (see text 
for further details on this choice). 



(a) (b) 




FIG. 3: (Color online): Purity [panel (a)] and logarithmic negativity 
[panel (b)] of the state of a two-site chain initially prepared in the 

maximally entangled state (l/\/2)(| } + | + +})i2 and affected 

by either thermal states of the oscillator environment (red dashed 
curve) or an entangled squeezed-vacuum state (solid dark one) for 
(HJ,hu},g) = (0.1, 0.01, 0.005)eV and r = 1.158, corresponding 
to a = 1/2. The revival at t — tt/uj associated with the squeezed- 
vacuum case is very pronounced. 
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FIG. 4: (Color online): Logarithmic negativity between the two sites 
of a chain affected by either thermal states of the oscillator environ- 
ment (green curve) or an entangled squeezed- vacuum state (red one) 
for r — 0.816 and a = 1/4. At t = n/uj, AT(t) has a revival, for 
the squeezed-vacuum case, which is absent from the thermal state- 
affected system. 



is valid for a squeezed-vacuum environment with 
A * = A^ q =n/(n + l). 

The additional oscillatory term proportional to A sq cos u>t 
appearing in Eq. (fl~8T > is responsible for the differences in the 
trend followed by the coherences corresponding to such two 
cases. The state purity is given by 

P(t) = \ + 2\(+ + \p s (t)\-~)\ 2 , (19) 

which is reported in Fig. [3] (a) for the two different types 
of environmental states considered in this study, highlighting 
in a rather striking way how strongly the entanglement be- 
tween the environmental oscillators affects the TLS's proper- 
ties. The purity revival at t — tt/lu is also reflected in the be- 
havior of the entan glem ent shared by the TLS's. We use log- 
arithmic negativity lfl3f1 as our entanglement measure. This is 
defined as Af(t) = log 2 ||pf2 T || with pf 2 T me partially trans- 
posed density matrix of the two TLS's and \ \A\\ = Tr yA^A 
the trace norm of any operator A. At t = tt/u>, Af(t) goes 
from zero (the system experiences a sudden death of entan- 
glement) to a rather large value. 

The trend illustrated so far has features of independence 
from the specific instance of initial state considered for the 
dimer. We have studied the evolution of Af(t) when the sys- 
tem is initialized in a separable state, finding qualitatively 
analogous behaviors of such figure of merit. For instance, 
Fig. |4] reports the results corresponding to a chain prepared 
as (l/\/2)(|— ) + |+))i ® \—)2 : entanglement builds up in the 
system (both the curves in Fig.|4]are such that Af(0) = 0) and 
never disappears (we have considered a relatively small cou- 
pling between environments and chain). However, while the 
case corresponding to a thermal-state environment has a peri- 
odicity of 2tt/uj, when a squeezing is present, a resurgence of 
entanglement occurs at half the thermal-state period. 

Such effects can be linked to an enhanced non-Markovian 
nature of the TLS dynamics arising when considering entan- 
gled environmental oscillators. The analysis of this connec- 
tion is the main aim of the following section. 



A. Characterizing the non-Markovian nature of the dynamics 

The quantification of the degree of non-Markovianity of 
adynamical evolution has been the focus of a significant re- 
cent research activity lfl6l - tl9ll . The main goal of such investi- 
gation has been the identification of appropriate tools for the 
characterization of the many facets of non-Markovianity and 
the proposal of operational ways to quantify it. More recently, 
such tools have found application in the characterization of 
the memory-keeping evolution of systems of physical rele- 
vance ll20l, l2lll . In this Subsection we move along such lines 
and employ one of the above-mentioned measures to deter- 
mine the features of the dynamics under study here and the 
differences between the two environmental configurations ad- 
dressed so far. 

Quantitatively, we take inspiration from the measure of 
non-Markovianity proposed by Rivasefa/. inRef. iTTvSl , which 
is based on the following approach: consider a bipartite en- 
tangled state such that only one component of the biparti- 
tion is affected by the dynamics we would like to charac- 
terize. Then, by computing the amount of entanglement be- 
tween the two parties at different instants of times within a 
selected interval [£rj,*max] we can detect non-markovianity 
by checking for a non-monotonic behavior of the quantum 
correlations. That is, for psa(0) = |$)(3>| with |<E>), some 
maximally entangled pure bipartite state, we define AE = 
E[psA(to)} - E[p SA (t max )} (where E denotes a legitimate 
entanglement monotone). We thus take 



r 


dE[p SA {t)\ 


10 


dt 



If the evolution of the system is Markovian, the derivative of 
E[p SA (t)] is negative andl^ = 0. 

While the measure in lfl7ll has been explicitly formulated 
for two qubits, so as to characterize the nature of a single-qubit 
channel, one can certainly think about extending the argument 
so as to address the case of a larger system. In particular, 
here we would like to be able to qualitatively understand the 
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0.5 1.0 1.5 2.0 

FIG. 5: (Color online): Witness of non-Markovianity for the thermal 
(red dashed line) and squeezed blue solid line) case for fuv = 0.01 
and r = 1.15 against the coupling strength a. 

reasons behind the occurrence of the resurgence peak when an 
entangled environment is at hand. 

Without embarking into the challenge of a rigorous formu- 
lation of the measure for non-Markovianity, we believe it is 
reasonable to consider the following situation: we take a gen- 
uinely tripartite entangled state of the two sites of the dimer 
at hand and an additional, ancillary TLS, which is exposed 
to a perfectly unitary dynamics (i.e. no environmental boson 
is interacting with the ancilla). More specifically, we con- 
sider the W state (1/V3)(| + - -) + 1- + -) + 1- - +))oia, 
which also ensures the presence of entanglement in any two- 
TLS reduction. While, as stated above, the ancilla a under- 
goes a unitary evolution, particles 1 and 2 are exposed to the 
effective quantum channel described so far. We then apply 
Eq. d20l) . where E is taken to be the logarithmic negativity 
between a and the dimer, to the cases where the environment 
is prepared in a squeezed-vacuum and a two-mode thermal 
state. The results are shown in Fig. [5] against the effective 
coupling rate a = g/hu> for £ max = 2tt/lj. While in the case 
of a thermal environment, I grows monotonically and satu- 
rates for a > 0.4, for a squeezed- vacuum state it peaks at 
about this value of a, signaling an increased non-Markovian 
behavior associated with such an environmental state. This 
magnified non-Markovianity is at the origin of the recurrence 
peaks exhibited in the analyses of the previous Sections. For 
lager coupling strengths a, each TLS -environment interaction 
becomes so strong with respect to the inter-TLS one that the 
system fragments into a series of local sub-systems, each be- 
ing basically oblivious to any non-local nature of the environ- 
mental state. 



III. CHAIN DYNAMICS UNDER A COMMON 
ENVIRONMENT 

We now move to the study of a configuration where the TLS 
chain is affected by a common environment. The expression 
for the system-environment coupling given in the first line of 
Eq. d2J is t nus modified into 

n a N 

itsR = g^2(a k + al)®J2 aZ j (21) 

k=l j=l 



where we have assumed an oscillator environment consisting 
of n a modes, each identically coupled to all the TLS of and 
TV-element chain. As remarked above, the symmetries in the 
present case are such that one can easily gather the exact an- 
alytical form of the evolved density matrix of the chain. For 
illustrative purposes, and to make a clear comparison to the 
case treated in previous Sections, we address a chain of two 
identical element, this time affected by a common environ- 
ment composed by two boson modes. In Ref. JUl, an anal- 
ogous problem has been approached by restricting the anal- 
ysis to the single-excitation subspace and using a numerical 
method for the analysis of non-Markovian dynamics. Here, 
we shall pursue an analytic approach that is made possible by 
the analysis conducted in Sec. Q] We start by considering an 
arbitrary pure state 1 77) = (a | — ) + at site 1 of the chain, 

while the second element is prepared in |— )2- The time evo- 
lution operator in the excitation-added displaced basis of the 
environmental modes is given by 

U(t) = 4 (/(t) ± | ± ±)(± ± I ® \<*t%)(ai%\ 

n,m 

+ [9(t)\ T ±>(T ± I + h(t)\ T ±)(± T |] ® \n,m){n,m\) 

(22) 




FIG. 6: (Color online) Purity against dimensionless interaction time 
ujt for thermal (red) and squeezed (green) environment states, with 
hj = 0.1, hw — 0.03, g — 0.01, a = 1/3 . In the two cases, 
squeezing parameter and (effective) temperature are r — 0.62 and 
T = 293 K. The dashed curves refer to the case of common environ- 
ments. 



(a) (b) 

jV(t) Af{t) 




FIG. 7: (Color online): Logarithmic Negativity for hj/g — 0.2 and 
a — 1/3, in the common- (dashed curve) and independent- (solid 
line) environment case. Panel (a) addresses the case of a thermal 
environment, while panel (b) is for a squeezed one. 
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with (p±(t) = e ±i ( 2e + ua8 ) t , i9(t) = cos(Jt) and n(t) = 
i sin( Jt). Notice that, if a = 0, then = |n). With this 



Very similar features between the individual-environment 
and the common-environment cases are found by studying the 
behavior of purity. Fig. [6] shows that, regardless of the en- 
vironmental configuration (whether made out of a common 
system or two individual oscillators), the squeezed-vacuum 
preparation is responsible for the revival of purity discussed 
above, a feature that is absent from the thermal case. 

It is interesting to notice that the central block of the density 
matrix in Eq. d23l does not contain any decoherence factor, 
thus being fully decoherence-free: entanglement and purity 
will be preserved in time, in such a subspace, which is associ- 
ated with a single excitation. Clearly, this marks a net differ- 
ence with the case of two independent environments, where 
no such special symmetry exists. Such a distinction is evi- 
dent from the inspection of Fig. [7] where we plot the loga- 
rithmic negativity in the common and individual-environment 
cases, for both the thermal and squeezed-vaccum preparation: 
the long-time behavior in the individual-environment arrange- 
ment clearly manifests a decrease of Af(t) that is absent from 
the common-environment configuration. 



IV. THREE-SITE CHAIN: THE ROLE OF THE 
SHARING-STRUCTURE OF ENVIRONMENTAL 
ENTANGLEMENT 

After having analyzed in the previous Sections the basic 
building-block for any complex structure of interacting TLS's, 
we move to the study of an open three-site chain (a trimer) af- 
fected by individual environments. While, qualitatively, the 
result gathered in Sec. [II] will be fully confirmed, such a sce- 
nario offers us an opportunity to investigate whether the shar- 
ing structure of the environmental entanglement plays any role 
in the open-system dynamics undergone by the TLS-chain. 

To start with, we have considered the three-mode squeezed 

state |sq 3 )=(l-A 2 )- 1 / 2 X;^ =0 A i |/,Z,Oi23- This state is 



expression at hand, the reduced density matrix of the chain is 



(23) 



(24) 



(25) 



such that, upon tracing out one of the modes, the remain- 
ing two are fully separable thermal states, thus witnessing the 
GHZ-like nature of such state. We have found negligible dif- 
ferences in the dynamics of the trimer arising from the use of 
such environmental state or the tensor product <E>j = iPth,j- as 
shown in Fig. [8] for instance, the purity function P(t) of the 
three-TLS state corresponding to such cases is very similar 
[as seen by inspecting the two bottom curves in Fig. [8). The 
insensitivity to squeezing is confirmed by looking at the bipar- 
tite entanglement between any pair of chain elements achieved 
after tracing out one TLS -oscillator subsystem. Fig. [9] illus- 
trates the behavior of the logarithmic negativity for such bi- 
partite states; the results are valid, qualitatively, regardless of 
the pair being taken and show that there is no revival peak at 
= tt/u in the case of a squeezed environment. Moreover, a 
similar behavior is observed when considering the tripartite 
negativity, which is a legitimate entanglement monotone for 
the tripartite case rf22"tl . 



P(t) 




FIG. 8: (Color online): Purity of three sites chain for the cases of 
thermal (green), generalized squeezed (red) and double split (blue) 
environment states. The nature of genuine tripartite entanglement of 
the generalized squeezed state is almost invisible, while some differ- 
ences in the time evolution can be pointed out for the split squeezed 
state. 



Pl2(t) = 



/0 

b 2 ti 2 {t) 
b 2 ${t)K{t) 



\ 

b 2 $(t)K*{t) ab*tp + (t)ti(t)i(t) 



b 2 n 2 (t) ab*<p+(t)K,(t)j(t) 



\0 a*b(p-(t)^(t)-y*(t) a*bip-(t)K*(t)j*(t) 
where we have introduced the decoherence factor 

7(t) = cxp[-2a 2 n(e l " t - i)( e -*"*+#" - 1)] 

for the thermal state case and 



•y(t) =cxp[— 2a 
for squeezed-vacuum one. 



2 (1 - e-' lut )(l - e"^*A 2 q ) + [1 - 2 cos(wt) + cos(2wf)]A sq 

AV 7 ! 1 
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The situation however changes if we modify the 
entanglement-sharing structure among the environmental os- 
cillators. To show this fact, we have taken a three-mode state 
obtained by superimposing at a beam-splitter one mode of a 
squeezed state |sq 2 ) 12 to an ancilla prepared in the vacuum 
state 1 0) 3 . More formally, we consider the state 

|sq 3 ) = (l! <g> B 23 )(9) |sq 2 ) 12 <g> |0} 3 (26) 

withB(8) = exp [(0/2) (q^as — a£a2)1 being the beam splitter 
operator of transmittivity cos 2 (6>/2) - F° r 9 = (9 = n), 
a two-mode squeezed-vacuum state of modes 1 and 2(1 and 
3) is retrieved, while for intermediate values of such param- 
eter, an entangled three-mode state is achieved having non- 
separable two-mode reduced states. In this second case, the 
environmental entanglement visibly affects the global proper- 
ties of the system in a way consistent with what we have seen 
earlier in Sec. [TTJ see Fig. [8] thus reinforcing the idea that the 
structure of quantum correlations shared by the environmen- 
tal oscillators plays a key role in the detailed non-Markovian 
dynamics of the TLS subsystem. 

V. CONCLUSIONS 

We have studied the full non-Markovian dynamics of a sys- 
tem of interacting TLS's arranged in a linear configuration and 
affected by a bosonic environment, considering both the indi- 
vidual and common-environment cases. By properly choos- 
ing the basis onto which each TLS-oscillator interaction is de- 
scribed, we have been able to push the analytic approach to 



such a problem up to the point of constructing a handy recipe 
for tracking any regime of interaction, from the weak to the 
strong-coupling regime with, in principle, an arbitrary num- 
ber of environmental modes per site. This has paved the way 
to the study of the effects that quantum entanglement shared 
by the elements of the environment has on the dynamics of 
the TLS system. Entanglement appears to enhance the non- 
Markovian character of the evolution, giving rise to specific 
features in relevant figures of merit of a TLS chain that are 
not present in the case of uncorrected environmental states. 
We have illustrated such effects by studying a dimer, which 
is the basic building-block for the construction of interesting 
interaction configurations that are usually considered for the 
description of biological complexes. Our investigation leaves 
open a series of questions related, in particular, to the scaling 
of the effects revealed here with the degree of connectivity, the 
depth and the geometry of the network of interacting TLS's. 
We plan to address such issues by extending and improving 
the implementability of our method. 
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